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Abstract. We study a toy model for "partially open" wave-mechanical system, like 
for instance a dielectric micro-cavity, in the semiclassical limit where ray dynamics 
is applicable. Our model is a quantized map on the 2-dimensional torus, with an 
additional damping at each time step, resulting in a subunitary propagator, or "damped 
quantum map". We obtain analogues of Weyl's laws for such maps in the semiclassical 
limit, and draw some more precise estimates when the classical dynamic is chaotic. 



1. Introduction 

A quantum billiard is a "closed quantum system", as it preserves probability. Math- 
ematically, this corresponds to the Laplace operator in Q with Dirichlet boundary condi- 
tions. In this case, the spectrum is discrete and the associated eigenfunctions are bound 
states. This system can be "opened" in various ways: among others, one possibility is to 
consider the situation where the refractive index takes two different values ?i-in/out inside 
and outside the billiard. This model can describe certain types of two-dimensional optical 
microresonators: after some approximations, the electromagnetic field satisfies the scalar 
Helmholtz equation (A -|- ^f^/^^J^ = inside and outside Q. For transverse magnetic 
polarization of the electromagnetic field, ^' and V^' are continuous across dQ, and the 
relation between kin/out and the energy E is expressed by k'fn/out ~ '^In/out^- ■^'^ ^^sA, 
case, the spectrum is purely absolutely continuous, and all bound states are replaced by 
metastable states: they correspond to complex generalized eigenvalues, called resonances^ 
which are the poles of the meromorphic continuation of the resolvent from the upper to 
the lower half plane. These quantum resonances play a physically significant role, as their 
imaginary part govern the decay in time of the metastable states. 

For such systems with a refractive index jump, the semiclassical (equivalently, the 
geometric optics) limit can be described as follows. A wavepacket travels along a single 
ray until in hits 50, then it generally splits between two rays, one reflected, the other 
refracted according to Snell's law. If the cavity O is convex, the refracted ray will escape 
to infinity, and we may concentrate to what happens inside: the wavepacket follows the 
same trajectory as in the case of the closed cavity, but it is damped at each bounce by 
a reflection factor depending on the incident angle. If we encode the classical dynamics 
inside the billiard by the bounce map, then the effect of that reflection factor is to damp 
the wavepacket after at each step (or bounce). This map does not preserve probability, 
it is a "weighted symplectic map", as studied in [NZj . 
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The toy model we will study below has the same characteristics as this bounce map, 
and has been primarily introduced in |KNS| to mimick the resonance spectra of dielectric 
microcavities. It consists in a symplectic smooth map k acting on a compact phase space 
(the 2-dimensional torus T^), plus a damping function on that phase space a G C°°(T^), 
with \a\ < 1. Throughout the paper we will not always assume precise dynamical property 
for the map k, althought the main result take a very specific form when k has the Anosov 
property. We will also assume that for each time n > 1, the fixed points of form a 
"thin" set (the precise condition is given in §2.21) . 

The corresponding quantum system will be contructed as follows: we first quantize 
K into a family of N x N unitary propagators, where the quantum dimension = 

{2-Kh)~^ =^ will be large, and the damping function a is quantized into an operator 
Op/j(a). All these quantities will be described in more detail in ^ To have a damping 
effect at the quantum level, we also need to assume that || Op/j(a)|| < 1 for all h <1. 

If we denote Uh{K) the unitary propagator obtained from a quantization of k, the 
damped quantum map then takes the form 

(1.1) Mh{a,K) = Oph{a)Uh{K). 

Apart from the ray-splitting situation described above, the above damped quantum map 
is also relevant as a toy model for the damped wave equation in a cavity or on a compact 
Riemannian manifold. Evolved through the damped wave equation, a wavepacket follows 
a geodesic at speed unity, and it is continuously damped along this trajectory [ALl |Sjo| . 
The above damped quantum map is a discrete time version of this type of evolution; it 
can be seen as a "stroboscopic" or "Poincare" map for such an evolution. To compare the 
spectrum of our quantum maps with the complex modes kn = uJn—i^ of a damped cavity, 
one should look at the modes contained in an interval |a;„ — /c| < vr around the frequency 
k PS : the distribution of the decay rates r„ : |a;„ — k\ < vr is expected to exhibit the 

same behavior as that of the decay rates {7!'*^ = — 21og|An'*^| , |Ara''^| G Spec(M/j(a, k))}, 
as we will see below. Some of the theorems we present here are analogues of theorems 
relative to the spectrum for the damped wave equation, proved in [ALl |Sjo| . Some of the 
latter theorems become trivial in the present framework, while the proofs of some others 
simplifies in the case of maps. Besides, the numerical diagonalization of finite matrices is 
simpler than that of wave operators. Also, it is easier to construct maps with pre-defined 
dynamical properties, than manifolds with pre-defined properties of the geodesic flow. 

We now come to our results concerning the maps (II. 1|) . Some of them - Theorems 11.21 
and 11.31 - have already been presented without proof in [NS] . 

In general, the matrix M/j(a, k) is not normal, and may not be diagonalizable. It is 
known that the spectrum of nonnormal matrices can be very sensitive to perturbations, 
leading to the more robust notion of pseudospectrum |ETj . We will show that, under the 
condition of nonvanishing damping factor \a\ > Omin > 0, the spectrum of Mh{a, k) is 
still rather constrained in the semiclassical limit: it resembles the spectrum of the unitary 
(undamped) map Uh{i^)- 
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Theorem 1.1. Let Mh{a,K) be the damped quantum map described above, where k is a 
smooth, symplectic map on and the damping factor a S C°°(T'^) satisfies 1 > |a| > 0. 

For each value of h = N~'^, N £ 'N, we denote by {X^'^^}j=i...h-^ the eigenvalues of 
Mh{a,n), counted with algebraic multiplicity. In the semiclassical limit h ^ 0, these 
eigenvalues are distributed as follows. Let us call 



n 



i=l 

and using the Birkhoff ergodic theorem, define 

Eloo(a) = ess inf lim a„ , ESoo(a) = ess sup lim . 

n— >oo n— >oo 

Then the spectrum semiclassically concentrates near an annulus delimited by the circles 
of radius Eloo(a) and ESoo(a)-' 

(1.2) V(5>0, lim /i#(l < j < /i^i : Eloo(a) -5 < lAf^l < ESoo(a) +^1 = 1. 

Suppose now that n is ergodic with respect to the Lebesgue measure ;U on T^, and 
denote 

(a) exp I / log |a|c?/i| the geometric mean of |a| on . 

In this case, the spectrum concentrates near the circle of radius (a) and the arguments 
of the eigenvalues become homogeneously distributed over §^ ~ [0, 1): 

Theorem 1.2. If k is ergodic with respect to the Lebesgue measure, 

{i) V(5 > 0, mil /i#|l < j < h-^ : \xf^\ - (a) < 6^ = I , 

h-^ are A'''^'* f 
(ii) yfGC\S'), lim/.^^/(^-^)=y^^/(t)dt. 

Note that (i) is an immediate consequence of Theorem 11.11 since the Birkhoff ergodic 
theorem states that Eloo(a) = ESoo(a) = (a) if n is ergodic with respect to /i. We remark 
that the ergodicity of k ensures that the spectrum of the unitary quantum maps Uh{K) 
become uniformly distributed as h ^ |BDB2l IMQK| . Actually, for this property to 
hold one only needs a weaker assumption, contained in Proposition 12.21 

If we now suppose that k has the Anosov property (which implies ergodicity), we can 
estimate more precisely the behavior of the spectrum as it concentrates around the circle 
of radius (a) in the semiclassical limit. 

Theorem 1.3. Suppose that k is Anosov. Then, for any e > 0, 

(1.3) lin, {l : I - (<.)| < (j^j^l = 1 ^ 
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This theorem is our main result, and is proven in ^ It rehes principally on the 
knowledge of the rate of convergence of the function a„ to (a) as n ^ oo, when k is 
Anosov. 

In general, for chaotic maps the radial distribution of the spectrum around (a) does not 
shrink to in the semiclassical limit, as numerical and analytical studies indicate [AL| . 
Toward this direction, we estimate the number of "large" eigenvalues, i.e. the subset of 
the spectrum that stay at a finite distance c > from (a), as /i ^ 0: we show that this 
number is bounded by h'^~^, where < < 1 can be seen as a "fractal" exponent and 
depend on both a and k. One more time, we use for this purpose information about the 
"probability" for the function to take values away from (a), as n becomes large. This 
involves large deviations properties for a„, which are usually expressed in terms of a rate 
function / > (see ^ depending on both a and k. For d > 0, considering the interval 
[log(a) + d,oo[, one has 

lim — log fi{x : log an{x) > d + log(a)} = —I{d) . 

n— >oo fi 

Our last result then the takes the form: 

Theorem 1.4. Suppose as before that k is Anosov, and choose c > 0. Define 



r = log(sup IIDkI^,'!!) , £c = log(l + c/(a)) , a_=rnin|a|. 

X 



T2 



For any constant C > and e > arbitrary small, denote C = C ±e, and set Ta k — 
(2r — 121oga_)~^. Then, if I denotes the rate function associated to a and k, we have 

< J < h-^ : \xf^\ > (a) + c} = ©(/I'^-^W) , 
where.. Ac) = ,i%f{^l . 

At this time, we do not know if the upper bound in Theorem II. 41 is optimal. The proof 
of the preceding result involves the use of an evolution time of order Ur ~ rlog/i~^, 
similar to an Ehrenfest time, up to which the quantum to classical correspondance - also 
known as Egorov theorem - is valid. The above bound is optimal in the sense that a 
particular choice of r = Tc makes minimal the bound we can obtain. It is given by 

def ^a,K 



a,K 



l + I{ic-)Ta. 

It is remarkable that in our setting, a small Ehrenfest time does not gives any relevant 
bound, but a large Ehrenfest time does not provide an optimal bound either, because the 
remainder terms in the Egorov theorem become too large. 

We also find interesting to note that in the context of the damped wave equation, a 
result equivalent to Theorems ll.ll and ll.2l was obtained by Sjostrand under the assump- 
tion that the geodesic flow is ergodic [SjoJ. But to our knowledge, no results similar to 
Theorem 1 1.31 are known in this framework. Concerning Theorem 11.41 a comparable result 
has been announced very recently in the case of the damped wave equation on manifolds 
of negative curvature [Ana], but working with flows on manifolds add technical compli- 
cations compared to our framework. Although, it could be appealing to compare the 
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nature of the upper bounds obtained in these two formahsms, and in both cases, it is 
still an open question to know if any lower bound could be determined for the number of 
eigenvalues larger than (o) + c, in the semiclassical limit. We also note that this "fractal 
Weyl law" is different from the fractal law for resonances presented in [SZ] . although the 
two systems share some similarities. 

In the whole paper, we discuss quantum maps on the 2-torus T^. The generalization to 
the 2n— dimensional torus, or any "reasonable" compact phase space does not present any 
new difficulty, provided a quantization can be constructed on it, in the spirit of [MOK| . 

In section [2l we introduce the general setting of quantum mechanics and pseudodif- 
ferential calculus on the torus. Theorems 11.11 and 11.21 together with some intermediate 
results are proved in section [Sj while the Anosov case is treated in section ID In section 
El we present numerical calculations of the spectrum of such maps to illustrate theorems 
11.11 ll.2l and ll.3l The observable a is chosen somewhat arbitrarily (with \a\ > 0), while k 
is a well-known perturbed cat map. 



We briefly recall the setting of quantum mechanics on the 2-torus. We refer to the 
literature for a more detailed presentation (HBI IBDBll [PEG] . 

2.1. The quantum torus. When the classical phase space is the torus =^ {x = 
{q,p) G (M/Z)^}, one can define a corresponding quantum space by imposing periodicity 
conditions in position and momentum on wave functions. When Planck's constant takes 
the discrete values h = {2ttN)^^, e N, these conditions yield a subspace of finite 
dimension which we will denote hy TCn- This space can be equipped with a "natural" 
hermitian scalar product. 

We begin by fixing the notations for the ^-Fourier transform on M, which maps position 
to momentum. Let S denote the Schwartz space of functions and S' its dual, i.e. the 
space of tempered distributions. The ^-Fourier transform of any ■0 G 5'(M) is defined as 



A wave function on the torus is a distribution periodic in both position and momentum: 



Such distributions can be nontrivial iS h = {2itN)^^ for some A^ G N, in which case they 
form a subspace TLn oi dimension A^. For simplicity we will take here 9i = 62 = 0. Then 
this space admits a "position" basis {\ej) : j £ Z/A^Z}, where 



2. Quantum mechanics on the torus 




ij{q + 1) = e^'^''i;{q), Fn^P{p + 1) = e^'^''Fn^{p) • 



(2.1) 




A "natural" hermitian product on TCn makes this basis orthonormal: 
(2.2) {ej\ek) = 5jk, j,k€Z/NZ, 

and we will denote by || • || the corresponding norm. 
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Let us now describe the quantization of observables on the torus. We start from 
pseudodifferential operators on L^(R) |GSl IDS| . To any / S 5(T*M) is associated its 
h—Weyl quantization, that is the operator acting on tp £ as: 

(2.3) /^^(g) ^ J fil±L^p)eiii-r)p^(r)drdp. 

This defines a continuous mapping from 5 to 5, hence from S' to S' by duahty. Further- 
more, it can be shown that the mapping f ^ can be extended to any / G C^(T*]R), 
the space of smooth functions with bounded derivatives, and the Calderon-Vaillancourt 
theorem shows that is also continuous on L^(M). 

A complex valued observable on the torus / G C°°(T^) can be identified with a biperi- 
odic function on (for all q,p, f{q+ l,p) = f{q,p+ 1) = f{q,p))- When h= l/27r7V, 
one can check that the operator maps the subspace Wat C 5'(R) to itself. In the 

def 

following, we will always adopt the notation h = 2nh, so that on the torus we have 
h = N~'^. This number h will play the role of a small parameter and remind us the stan- 
dard ^-pseudodifferential calculus in r*M. We will then write Opi^{f) for the restriction 
of on Ti-N, which will be the quantization of / on the torus. It is a iV x iV matrix in 
the basis ([O) . 

The operator Op/j(/) inherits some properties from f^. We will list the ones which will 
be useful to us. Op^(/)^ = Op^(/*), so if / takes real values, Op/^(/) is self-adjoint. The 
function f = fh may also depend on h, and to keep on the torus the main features of the 
standard pseudodifferential calculus in T*R, these functions - or symbols - must belong to 
particular classes. On the torus, these different classes are defined exactly as for symbols 
in C^{T*W). For a sequence of functions {h)He]o,i], fh = f{x,K) £ C^{T*R) x]0,l], 
we will say that fh G 5^(1), with 6 g]0, m G R if h^ffi is uniformly bounded with 
respect to H and for any multi-index a = (ni, 7^2) G of length \a\ = ni + 71-2, we have : 

In the latter equation, || • \\co denotes the sup-norm on and 9" stands for ^TifJ^- 
On the torus, one simply has 27r?i = h = N^'^ for some G N, and C^(T*M) is replaced 
by C°°(T^). Let us denote these symbol classes. We have the following inequality 
of norms, useful to carry properties of ^—pseudodifferential operators on R to the torus 
[BDBlj : 

(2.4) V/gC°^(t2), ||Op,(/)||<||/^k.^i.. 

Note that this property remains valid if / = //j depends on h, with fn, G 5'^(1). The 
continuity states that if fti G ^^(l), then is a bounded operator (with ^—uniform 
bound) from L^{R) to L^{R). Since Ophif) is the restriction of on Hn, Eq. M 
implies the existence of a constant C independent of h such that for / G 5"^ and h g]0, 1], 



Op,(/)||<C. 
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The symbol calculus on the Weyl operators f '^ easily extends to their restrictions on 
Ti-N- We have, for symbols / and g in 5^ the composition rule: 

0Phif)0pdg) = 0pM9), 
where f^hd is defined as usual: for X = (x,^) G we have 

(2.5) (/MW = j^J{X + Z)g{X + Y)e't-^^^^)dYdZ 

where a denotes the usual symplectic form. Another useful expression for f^hQ is given 
by 

(2.6) (/M(^) = e^'^(^-^-)(/(X)<7(y))|x=Y 
where Dx = {Dx, D^) = {jdx, jd^). This representation is particulary useful for /i— expansions 

T and g E S^, then fho S S'p 



of f\lhg- If / e and g € S?, then /Jj/^g e SV''^"', and at first order we have : 



(2.7) II Op,(/) Op,{g) - Op,{fg)\\ < Cf,gh'-''-^-- . 

def 

The Weyl operators on Hn, Tm,n = OPhi^mn), with 
allow us to represent Op f^{f): 

(2.8) Op;j(/) = ^ fm,nTm,n, where fm,n = / f dji . 
From the trace identities 

^2 9^ _ ^{-ly•'^"' h'^ if/i,zy = mod/i-i 

^''^ 1 otherwise, 

one easily shows that for any f £ Sg we have 

(2.10) /iTV(Op^(/)) = /o,o + 0(/i~)= / fdfi + Oih^). 



Let a G ^q. We will write a_ = minx2 |a|, and a+ = maxx2 |a|. The next proposition 
adapts the sharp Carding inequality to the torus setting. 

Proposition 2.1. Let a £ he a real, positive symbol, with Rana = [a_,a+]. There 
exist a constant C > such that, for small enough h and any normalized state u £ 7iN-' 

a_ - Ch^-^^ < {u, Oph(,a)u) < a+ + Ch^-^^ . 

Proof. We first sketch the proof of the sharp Garding inequality in the case of pseudodif- 
ferential operators on T*]R with real symbol a G 5'^(1). For the lower bound, we suppose 
without loss of generality that a_ = 0. 

For X = {x,0 e consider T{X) = T{x,0 = 2e"H^. Using Eq. ([231) for 
symbols on T*]R, a straightforward calculation involving Gaussian integrals shows that 



EMMANUEL SCHENCK 



rft/iF = r, hence is an orthogonal projector, and thus a positive operator. Define the 
symbol 



a*r(X)=^ — -/ a{X + Y)2e~^ dY . 



To connect (a*r)^ with a^, we make use of the Taylor formula at point X in the above 
definition. Because of the parity of F and the fact that J^2 T{X)dX = 2Trh, we have 

a*r{X) = a{X) + ^ l'l'{l-e)a"{X + eY)Y^ e'^ dVde 



aiX) + i - 0)ha'{X + ezVh)Z^ e"^" dZd9 



r(X) 

To evaluate ||?^^||l2_>j^2, we rescale the variable X ^ X = X/^/h, and call f{X) = 

r{VhX). If M G we denote u{x) = Hiu{x). This transformation is unitary : \\u\\l'2 = 
1 1 It 1 1 £,2. Now, a simple change of variables using (|2.3I1 shows that Hr^uH^.a = ||^J"'u||2,2 
where denote the h = 1 quantization of the symbol f. Because of the term ha" = 
0{h^~^^) appearing in the definition of r, for any multi-index \a\ we have 

Since the continuity theorem applied to rf yields to a bound that involves a finite 
number of derivatives of f, we get H^^f ||l2-»l2 = 0{h^~'^^). Here and below, by / < 5 we 
will mean that |/| < C\g\ for some C > 0. Hence, 

and we conclude by ||?^^||l2_^^2 = 0{h}^'^^) . 
Now, from the definition of a^T, we also have 



(«*r));; = — y a{Y){T{. -YmdY . 

But as we noticed above, (r(- — 1"))^ > 0, and then (a-kV)'^ is positive definite. Using 
the fact that a -kV = a + r and ||r^||i2_^/,2 = 0{}t}~'^^), this implies the existence of a 
constant c > such that 

(2.11) (n, a^u) > -c/i^-^'^ . 

The upper bound is obtained similarly, assuming a+ = and considering the symbol 
-a > 0. 

It is now a straightforward calculation to transpose these properties on the torus by 
making use of Eq. (12. 4p . □ 



WEYL LAWS FOR PARTIALLY OPEN QUANTUM MAPS 



9 



2.2. Quantum dynamics. The classical dynamics will simply be given by a smooth 
symplectic diffeomorphism k : ^ T^. Since we are mainly interested in the case 
of chaotic dynamics, we will sometimes make the hypothesis that k is Anosov, hence 
ergodic with respect to the Lebesgue measure ^. From this ergodicity we draw the 
following consequence on the set of periodic points. We recall that a set has Minkowski 
content zero if, for any e > 0, it can be covered by equiradial Euclidean balls of total 
measure less than e. 

Proposition 2.2. Assume the diffeomorphism k is ergodic w.r.to the Lebesgue measure. 
Then, for any n > 1, the fixed points of k" form a set of Minkowski content zero. 

Proof. Let us first check that, for any n / 0, the set of n-periodic points Fix(K'^) has 
Lebesgue measure zero. Indeed, this set is K-invariant, so by ergodicity it has measure 
or 1. In the latter case, the map n would be n-periodic on a set of full measure, and 
therefore not ergodic. 

Let us now fix n 7^ 0. Since is continuous, for any e > the set 

=^ {x e T^ dist(x, < e} is closed in T^. 

Since C F^i if e < e', for any Borel measure u on we have 

v{Fq) = lim v{F,) . 

Since Fq = Fix(K'^) has zero Lebesgue measure, it is of Minkowski content zero. □ 

We will not study in detail the possible quantization recipes of the symplectic map 
K (see e.g. [KMRl iDEGl IZel] for discussions on this question), but assume that some 
quantization can be constructed. In dimension d = 1, the map k can be decomposed into 
the product of three maps L, tv and where L G SL(2, Z) is a linear automorphism of 
the torus, ty, is the translation of vector v and (pi is a time 1 hamiltonian flow |[KMR| . One 
quantizes the map k = Lot^ocpi by quantizing separately L, tv and (pi into Uh{L), Uh{t^) 
and Uh{(pi) and setting Uh{n) = Uh{L) Uh{tv) Uh{(pi) ■ We are mainly interested in the 
Egorov property, or "quantum to classical correspondence principle" of such maps, which 
is expressed for / G 5^ by 

\\Uh{K)-' Opj,{f )Uh{K) - Op^^if o k)\\ = Oh{l) , where Oh{l) ^ . 
The following lemma makes the defect term Oh{l) in the preceding equation more precise. 
Lemma 2.3. Let f G S^. There is a constant Cj^i^ such that 
(2.12) WUhif^r' Opj,{f)Uh{K) - OpJ/ o k)|| < Cf,^h'-^' . 

Proof. Since it is well known that for linear maps L, one has 

Uh{Ly^ Ovh{f)Uh{L) = Ophif ° L) , 

we will consider only the quantizations of iv and (/>i. The map t^ is quantized by a 
quantum translation operator of vector v/^, which is at distance |v — v/j| = 0(/i): 

UhiUr^ Ovh{f)Uh{t^) = Ovhif o tvj . 
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Consequently, we need to estimate || Opf^{f o — / o tv^)!!- For this purpose, we first 
evaluate ||/ o — / ° ||co- Denote v = {v'^^v^) and v/j = [vf^^v^^). A Taylor expansion 
shows that 

\\f{q + vl + {v'^ - vl),p + vl + K - vD) - fiq + vlp + vl)\\co = 0{h^-') 
and then, 

||/otv-/oivJlc0=O(/i'~'). 

Hence, h^^^{f o — / o t^^) g Sg. Using Proposition 12.11 we conclude by 

\\0pf,{fou-fou^)\\=0{h'-'). 

Let us denote hy H Sq the Hamiltonian, Xh the associated Hamiltonian vector field, 
and (pt = exp (iXn) the classical Hamiltonian flow at time t. For simplicity, we denote 

the quantum propagator at time t by '= exp(— ^ Opf^{H)), and write 0^ / *== / o (p^. 
In particular, note that = Uh{(pi)- From the equations: 

irt+sf\s=o = {H,rtf} 

j^U'+^ Ovy,{4>UJ)U-^'+%=^ = Difft(F, f)u-' 
with Ty\UHJ) = Ov^,{{HAU}) - k[^Vh{H),OpMtf)l we get 

OpMf) =U-'Op^{f)U'+ ['w-'BiSsiH,f)U'-'ds. 

Jo 

A straightforward application of (|2.6p yields to 

^[Op^{H),OpM.f)] = Op </.?/}) + ^OnAh'-'') , 

where O-j^j^ {q) denotes an operator in Ti.N whose norm is of order q. If we use the unitarity 
of U^, we thus obtain : 

\\U-' Op,{f)U' - OpMf)\\ = o{h'-'') . 
Adding all these estimates, we end up with 

WUhiK)-' Opj,{f)UhiK) - Opj^if ok)\\< Cf,^-^' . 

□ 

Taking into account the damping, our damped quantum map is given by the matrix 
Mfi{a, k) in The damping factor a G 5q is chosen such that, for h small enough, 

11 Op;^(a)|| < 1 and a_ = m.mj2 \a\ > 0. From Proposition 12.11 this implies that Op;,(a), 
and thus Mh{a,K), are invertible, with inverses uniformly bounded with respect to h: 

(2.13) \\Mh{a,K)\\ = a+ + 0{h), \\Mh{a, k)-'^\\ = aZ^ + 0{h) . 

As explained in the introduction, Mh{a, k) is not a normal operator, and it may not be 
diagonalizable. Nevertheless, we may write its spectrum as 

Spec(M,(a, k)) = {aS"\ A^"\ . . . , Ai^'A} , 
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where each eigenvalue is counted according to its algebraic multiplicity, and eigenvalues 
are ordered by decreasing modulus (in the following we will sometimes omit the ^'^^ 
supersripts). The bounds (|2.13p trivially imply 

(2.14) a+ + Ch> lAf^l > \xi^^\ > ... > |Af_\| >a.-Ch, 

for some constant C > 0. Since we assumed a_ > 0, the spectrum is localized in an 
annulus for h small enough. The above bounds are similar with the ones obtained for 
the damped wave equation [ALi Eq.(2-2)]. 

For later use, let us now recall the Weyl inequalities |Konj . which relate the eigenvalues 
of an operator A to its singular values (that is, the eigenvalues of V A^A): 

Proposition 2.4 (Weyl's inequalities). LetTL he a Hilbert space and A G C{Ti.) a compact 
operator. Denote respectively ai,a2,... (respectively Pi,P2,---) Us eigenvalues (resp. 
singular values) ordered by decreasing moduli and counted with algebraic multiplicities. 
Then, VA; < dimTi, we have : 

k k k k 

(2.15) ni«^i^n/^^' 

1=1 i=l i=l i=l 

We immediately deduce from this the following 

(n) (n) 

Corollary 2.5. Fix n > 1. Let Ai, A2, ■ ■ ■ be the eigenvalues of A and s\ , • • • the 

eigenvalues of An = ^a/ A^^A^, ordered as above. Then for any k < diuiTi: 

x:iogiA.i<x:iog.i"\ fiiA.i<n.i"). 

i=l i=l 1=1 i=l 

2.3. A first-order functional calculus on C{TLn)- In iSwe will need to analyze the 
operators 

{yA4(a,K)t" A4(a,/€)", n G N} . 

We will show that they are (self-adjoint) pseudodifferential operators (i.e. "quantum 
observables") on TLn, and then draw some estimates on their spectra in the semiclassical 
limit via counting functions and trace methods. We will use for this a functional calculus 
for operators in C{TLn), obtained from a Cauchy formula, via the method of almost 
analytic extensions. 

Let a G be a real symbol. In order to localize the spectrum of Op/j(a) over a set 
depending on /i, we will make use of compactly supported functions /«, G Co"(M) which 
can have variations of order 1 over distances of order w{h), for some continuous function 

w > 0, satisfying w{h) Q. 

The functions fw have derivatives growing as /i — > 0: for any m G N, we will assume 
that 

(2.16) \\d"'fn.\\co<Cmwih)-"'. 
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We will call such functions w{h)— admissible. Our main goal in this section consist in 
defining the operators fw{Opfi{a)) and characterize them as pseudodifferential operators 
on the torus. 

We first construct an almost analytic extension fw G C^(C) satisfying 



\\dU\co < c^\Qzrw{hy 



-m-2 



Vm > 



(2.17) 

(2-18) = fyj , 

where d stands for ^ = \[-§^ + i^)- For this purpose, we follow closely [DS] : one put 
X G Co°(I^) be equal to 1 in a neighborhood of 0, V' G C'(5"(]R) equal to 1 in a neighborhood 
of Supp(/^) and define 

~Ux + iy) = ^^^^ / eH^+'^)«x(y6/U0de. 

We will now check that fw is a function that satisfies Eqs. (I2.17P and (I2.18p . Notice first 
that fw is compactly supported in C, and that Eq. (|2.18ll follows clearly from the Fourier 
inversion formula. To derive (|2.17p . remark that 



\ilj{x)x{y) 



Air 



-y 



x'jyQ . 
{yO"" 



+ 



'>P'{x)x{y) +i'>P{x)x'{y) f 



47r 



xiyC)fwix)dxd^ 



I + II. 



Note that if y = 0, I = II = because of the properties of x and tp, so (|2.17p is satisfied 
with Cm = 0. We now suppose y 7^ 0. Since is compactly supported, we can integrate 
by parts and using (I2.16p . we obtain, setting t = y^ : 



x'jyQ . 
(yCY 



< 



x'{t) 



(dT^'Ux)) 



-j-m+l 



dtdx < Cmw{h) 



-m-2 



In the last inequality, we used the fact that and t 1— > t ™ ^x'{t) are compactly 
supported. Hence, this shows that |I| < Cm|y|™'u;(/i)~'"~^. To treat the term II, denote 

Fi^^^y) = ^'^^)^^y) + -^^^^)^^y) and G{x,x^y) = {:. + D,fDn- ^^^^^ ^ 



47r 



'x — X + iy 



Since y / 0, we can rewrite II to get 

jix-x+iy)^ x'jyQy 
i^ + ^nyO' 



II = iF{x, y) J y™ e 
As above, we set t = y^. This gives 



;G{x, X, y)dxd^ . 



-G(x,x,y) 



dxd£, < 



m+2 



x'{t) 



{t + iy)2f™ 



G(x,x,y) 



dtdx 



< 



y"'^'02Gix,x,y) 



dtdx . 



Let us distinguish two cases. 
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• If X ^ SuppV', then x G Supp/^i =^ x — x 7^ 0, from which we deduce that 
J\G{x,x,y)\dx<Cmw{h)-^-^. 

• If X E Supp-fA, then F{x,y) = i(47r)~^x'(y)- If x'iu) = 0; we get II = 0. Other- 
wise, we have bi < \y\ < 62 for some fixed constants 61,62 > depending on x- 
In this case, we get again / \G{x,x,y)\dx < Cmw{h)~"^~'^ . 

Grouping the results, we see that |II| < Cm\y\^'^'^w{h)^"^^'^ < Crn\y\"^w{h)~'^~'^ , and 
since |I| < Cmlyrwih)-"^-^ , it follows that Eq. (121711 holds. 

Considering a function f^j as above, we now characterize fw{OPhi^)) a pseudodif- 
ferential operator on the torus. We begin by two lemmas concerning resolvent estimates. 
For z ^ SpecOp;j(a), we denote 'R-z{a) = (Op/j(2; — a))~^ the resolvent of Op;j(a) at point 
z. 

Lemma 2.6. Let a £ Sg be a real symbol, and Q C C a bounded domain such that 
Supp(/«,) C Jl. Take z £ fl and suppose that \Qz\ > /i^ for some e €]0, 1] such that 
6 + e < Then, 



z — a 



Proof. The hypothesis \Qz\ > implies immediately that 

(2.19) r^-^<h~'. 

\z — a\ 

To control the derivatives of (2; — a) we will make use of the Faa di Bruno formula 
[Comj . For n > 2, Let a = {ai,...,an) G N" be a multi-index, and IT be the set of 
partitions of the ensemble {1, \a\}. For vr G 11, we write vr = {Bi, ...Br}, where Bi is 
some subset of {1, \a\} and can then be seen as a multi-index. Here |a| > r > 1, and 
we denote |7r| = r. For two smooth functions g : M" ^ M and / : M M such that fog 
is well defined, one has 

(2-20) 9Vog = Y:^d\-\f{9)Yl . 

We now take n = 2, and q G N^. Using this formula for 5" and recalling that a £ Sg, 
we get for each partition vr G 11 a sum of terms which can be written : 

^ TT 9^(z - a) < . 

(z - a V^+i ii ^ ^ ~ 
^ ' Ben 

Since we have r < |a|, this concludes the proof. □ 

The preceding lemma allows us to obtain now a useful resolvent estimate : 

Lemma 2.7. Choose e < Suppose as above that a £ Sg is real and z £ with 

> /i^ Then, 

7^,(a) = 0p^ (j^) +^^(^)' 
where Rh{z) £ C{Hn) satisfies \\Rh{z)\\ = C'(/i^~^('^+^'^)) , uniformly in z. 
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Proof. We will denote O-H^^zig) an operator which depends continuously on z and whose 
norm in TCn is of order q. By the preceding lemma and the symbolic calculus (12. 7p . for 

\Qz\ > h'^ we can write : 

Ophiz - a) (t") = - OPhir-z) 
with Ophirz) = C'hjv,^(^^"^'^"^^)- For ^ small enough, the right hand side is invertible : 

(2.21) (Id-Op^(r,))-i =Id+0„^,,(/ii-2^-3-). 

We now remark that the Garding inequality implies Op/j (^jhi) ~ ^'Hn,z{^^^) ■ Since 
we have obviously 

TZzia) = Op,, (j^) (Id - Op;,(r,))-i , 



we obtain using Eq. (I2.2ip 



l-2(5-4£\ 



7^,(a) = Op;, [t^j + ^-hnA^ 

It is now straightforward to check that all the remainder terms are uniform with respect 
to since we have z G and \'isz\ > . □ 

We can now formulate a first order functional calculus for a € S^. 

Proposition 2.8 (Functional calculus). Let a £ real, and a w— admissible function 
with w{h)~^ < h~'^ such that 

-I 2^ 

(2.22) < 7/ < 



6 

Then, for any e > such that rj < e < — g— , we have : 

/^(Opja)) = Op^(/^(a)) + OnAh^'^''^') 

and fw{a) G 

Proof. Let us write the Lebesgue measure dxdy = Prom the operator theory point 

of view, we already know that for any bounded self-adjoint operator A 

fUA) = ^ [ dfUz, z){z - A)-^ dzAdz, 

see for example [DSj for a proof. 

Because of the condition expressed by Eq. (|2.22|) . it is possible to choose e > such 
that rj <£ < We now divide the complex plane into two subsets depending on e: 

C = U where f^i = {z G C : < hf} and 1^2 = C \ fii. We first treat the case 
z e r^i. We recall that for z ^ R, one has 

(2-23) ||^.(«)ll < ^ and < . 

\xsz\ \z — a 
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Using Eqs. (|2.17D and (I2.23p . we obtain for any m G N 



dfw{z, z)lZz{a) dz A dz 



< 



1 



Cm\Qzrw{h)-'^-^dz A dz 



Since h'^w{h) ^ ^ ^ h^o^ taking m sufficiently large we obtain 



dfrnTlzia) dz A dz 



We now consider z ^VL2. Using Lemma [277l we obtain 

dz A dz 



dU{z,z)nz{a) 
Q.2 2i7r 



dfw{z,z) 



^2 



1 \ dz A dz 



(2.24) 



+ / dfw{z,z)Rh{z, 
Jn2 2i7r 

I + II. 



z — a J 2i7r 

dz A dz 



To rewrite the first term I in the right hand side, we use the Fourier representation (12.8(1 : 
I = ^ / dfi, ^ ff,^uC,,,u(.z)dz A dz 



Op, 



T^,u / dxdye^'^ 



(px-uy) _ 



dfwiz,z) 



2ivr z - a{x, y) 



dz A dz 



2i7r 



dfwiz,z] 



dz A dz 



^2 



z — a 



where we used the uniform convergence of the Fourier series (12.8(1 , the Pubini Theorem 
and the linearity of the quantization / ^ Opf^{f). Note that 



dfw{z,z) 



dz A dz 



dfw{z,z] 
c z-a 



dz A dz 



dfw{z,z] 



dz Adz . 



IQ2 ^ - a 
Now, Eqs. (12:2311 and (l2T7ll yields to 

f ?M^dzAdz = 0{h^), 
Jn^ z-a 

and since the Cauchy formula for C°° functions implies that 

1 f dfw{z,z) 



2i7r r z — a 



-dz Adz = fw\R{a) = fu,{a) 



we can simply write 
(2.25) I = Op;, 



2i7r 



dfw{z,z) 



n2 



z — a 



dzAdz] =OpJ/^(a)) + 0„^(/i~). 



< / \\Rh{z)\\\dU{z,z)\dzAdz 
Jn2 
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For the term II in Eq. (I2.24p . we first remark that dfw is compactly supported. Then, 
we use Eq. (|2.17ll and Lemma [2771 to write : 

dfw{z, z)Rh{z)dz A dz 

(2.26) < /ll-25-4.;^-2„ < ^1-25-6. 

Hence, combining Eqs. ^a?M . ^a?M and ^0^ . 

1 
2i^ 

as was to be shown. The last assertion is a direct application of Eqs. (I2.20p and (I2.16p . □ 



dU{z,z)nzia)dz Adz = Op ^ifw (a)) + On ^ih 



l-25-6£N 



Below, we will frequently make use of the functional calculus for some perturbed 
operators. Our main tool for this purpose is stated as follows: 

Corollary 2.9 (Functional calculus - perturbations). Let a ^ he a real symbol and 
fyj he a w— admissible function with w{h)^^ < and < i] < Consider E 

C{'Hn), with the properties : 

(i) Opfi{a) + Bfi is self-adjoint, 

(ii) \\Bh\\ = 0{h'^)for some u > Arj > . 

Then, at first order, the functional calculus is still valid : for any e > such that 

1-26 V 
77 < e < — - — and rj < e < — , 
6 4 

we have : 

/^(Op,(a) + Bh) = Op,(/^(a)) + 0«^(/7--(i-2^-6='--4=)) 
and fw{a) G S^^_^^. 

Proof. We must find the resolvent of Opj^{a) + B^. From {ii) we can choose e > such 
that rj < £ < and rj < e < j. Note that this implies v — ie > 0. As in Lemma 
12.71 we choose a compact domain Q C C with z G and split O = Qi U O2 as above. 
Suppose first that z £ 0,2, i-e. |3'(-z)| > ■ Then, using the condition (ii) we get 

Id-(0p;,(z-a)-i?;,)0p^(^) = 0p;,(r,) + i?;,0p/,(- ^ ^ 



z — a z — a 



The norm of the right hand side is of order 25 3£,u e) uniformly for z £ hence 

for h small enough, we can use the same method employed in the Lemma \2J\ to get : 



{z - Ophia) - Bh)-' = Op^{^)[ld-Opf,{r,)-BhOphi 



1 . / „ , . „ , 1 ^ 



z — a \ z — a' 



Op^,{^)+R, 
z — a 
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where \\Rh\\ ^ j^nim{i-2S~4e,u-2E) uniformly in z. The next steps are now exactly the 
same as above, and we end up with 

where now \\Rh{z)\\ < /i™n(i-25-4e,,.-2e) ^ gjj^pg P/^lbo < h-^"), this concludes the 
proof. □ 

As a direct application of the last two proposition, let us show an analogue of the 
spectral Weyl law on the torus. 

Proposition 2.10. Let a £ S^, and Bh G C{nN) as in Corollary \2.9[ Choose e: > 
arbitrary small but fixed, and call Ah = Op/j(a) + Bh Let Ei,E2 be positive numbers, 

I =^ [Ei,E2] and I± =^ [Ei :fe,E2±e]. Then, 

(2.27) / li_{a)dfi + Oh{l) <h#{XG Spec Ah nl} < [ l/+(o) + o;,(l) . 

Proof. Define a smooth function such that for some C > 0, = 1 if x G / and 

X^{x) = if X ^ /+. Define as well x~ = outside / and x~ = 1 on /_. Denote 
Dh = h#{X G Spec Ah D I}. Then, 

hTr{x-{Ah)) <Dh< hTT{x+{Ah)) . 

By the Corollary EH /iTr(x+(Op;,(a) + Bh)) = fj2 X^{a)dix + 0{h'^) for some a > 0, 
and /iTr(x~(Op^(a) + Bh)) = fj2 X~{o)diJL + 0{h^) as well. But obviously. 



This yields to 



/ ^I.{a)d^i< I x'{a) and [ x^{a)< [ 

Jj2 Jj2 Jy2 Jj2 

[ li_ia)df^ + Ohil)<Dh< [ ti4a)+Ohil). 
Jt2 7t2 



□ 



3. Eigenvalues density 



3.1. The operator 5„. Consider our damping function a G 5q, with Ran(|a|) = [a_,a+], 
a_ > 0, a+ < 1. To simplify the following analysis, we will suppose without any loss of 
generality that a+ = 1. As mentioned above, to study the radial distribution of Mh{a, k) 
it will be useful to first consider the sequence of operators 



(3.1) Ma) = Mhia, K)t" Mh{a, n > 1 . 

Let us show that for /i — > and n > 1 possibly depending on /i, these operators can 
be rewritten into a more simple form, involving n— time evolutions of the observable a 
by the map k. Using the composition of operators (|2.7|) and the Egorov property (|2.12|) . 
we will show that the quantum to classical correspondence is valid up to times of order 
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log/i^^, as it is usually expected. In what follows, for any constant C, we will make 
use of the notation = C ± e with e > arbitrary small but fixed as /i ^ 0. We 
allow the value of e to change from equation to equation, hence denotes any constant 
arbitrary close to C, being larger than C and C~ smaller: e will then be chosen small 
enough so that the equations where appear are satisfied. We also recall the following 
definitions, already introduced in Theorem ll.4t 

(3-2) r = log(sup ||L>kU||) and Ta,^ = — — -i- . 

X zi — 12. log a_ 

Proposition 3.1. LetT > be a constant such thatr < Ta^n- If E{x) denotes the integer 
part of X, define Uj- = E{Tlogh~^). Then, for n < n^, the operators 

iSn{a) = log5„(a) 
are well defined if and only ifKeT{Mh{a, k)) = 0. Furthermore, if we set 

def j~r I j I i „ def i i i i\ 
Un = IMOO'^I", (-CLn = — >, lOg | a o K | , 

1=1 ^ 1=1 

we have for n < n,-; 

(3.3) iSnia) = Opf,{ian) + Ouj, {h^' ) 

(3.4) Sn{a) = Op^(a„) + (h'^' ) 
where a = 1 — r/Ta^K > 0. 

Proof. Let us underline the main steps we will encounter below. Writing first 

n 

def "TT I 1 2 

CLn = I I ° «^ I 

i=l 

for n G N, we show that for r < a„ belongs to a symbol class 5^+ with 6^ < 1/2. 
Then, we show by using the symbolic calculus (12. 7p and the Egorov property (I2.12P that 
'Snia) = Opi^{an) + 0{h'^ ) for some > 0. Finally, by bounding the spectrum of 5„(a), 
we will complete the proof of the proposition by using the functional calculus to define 
and compute both iSn{a) and 5„(a). 

Lemma 3.2. Let a £ Sq be a symbol on the torus. Then, for any multi-index a G N^, 
there exists Ca,a,K > such that 

(3.5) Vn> 1, ||a"(aoK")||co < C,,a,.e"l°|r . 
Hence for any r > such that t < we have uniformly for n < n,-.- 

(3.6) aoK"G5^ 6 = Tr<^. 
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Proof. The behavior expressed by Eq. (|3.5|) is well known for flows |BR| . and we refer to 
[FNWj . Lemma 1 for a detailed proof in the case of applications. The second part of the 
lemma follows easily from (13. 5p . □ 



Lemma 3.3. If t < we have : 

yn < n-r , a„ G 5^+ with 5^ = tT^ < - . 

Proof. Since a+ = 1 and < a_ < 1, is uniformly bounded from above with respect 
to h. It is then enough to show that for every multi index a G and n < Ur, one has 



i=l 



Set by convention d^f = f. Applying the Leibniz rule, we can write 



ao K 



i=l 



< n 



\a\ 



sup 



CO 



«i + ...+a„=a 



[aoK )\\co 



Let us look at a typical term in the product appearing in the right hand side. Since at 
most I a I indices a,- in are non zero and la o < 1, 



i=l 



i=l 



< ( sup C/3,a,.)l"le"l-|r1^'ir„,,,,e"l"|r . 
\P\<H 



Finally, we simply get 



i=l 



< JC p|Q!|{rn+logn) ^ 7^ p|o|r'n 



CO 



for some T' > T. If we choose T' = r+ and r+ — F small enough such that tF+ < i, 
the last equation will be true only for n > uq, with no fixed independent of h. Since for 
n < no, we obviously have a„ G 5*0, we finally conclude that if tF+ < ^, a„ G 5^+ with 
= tT~^ < |, uniformly for n < n^-. □ 

We can now rewrite more explicitly Eq. (|3.ip . 

Lemma 3.4. Choose r > small enough such that 5 = tT < ^, and take as before 
n <nr. Then, we have 

(3.7) 5„(a) = Op;,(a„) + 

where u = 1 — 26. 

Proof. The Lemmas 13.21 and 1 3.31 tell us that uniformly for n < E{t log h^^), the symbols 
ook" and belong to the class 5^+ ii6 = tT < ^. This allows to write (using U = Uh{K) 
for simplicity): 
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(3.8) =U^OpM^)U + OnAh'-^'^) 

(3.9) =Op;,(|aoK|2)+0„^(/ii-25+) 

where the symbohc calculus (|2.7p has been used to get (|3.8p and Eq. (|2.12l) to deduce 
Eq. (13. 9p . The remainder R\ have a norm of order h^~'^^~^ in Hn- This calculation can 
be iterated : suppose that the preceding step gave 



k 

k 

i 1 

i=l 



with i?f = O-Hj^ih^-'^^^) for 1 < i < A;. Then, with the same arguments as those we used 
for the first step, we can find R^^\ = 0-)^j^{h^~'^^^) such that 

Ophia) Opjafc) Oph{a)U = Op/,(afc+i) + R^ . 

If we define now = Opf,{a)R^ Opf,ia)U = Onj,{h^'^^^ ), we get 

k 

U^Opj,{a)Sk{a)Opj,{a)U = Opj,{a) Opj,{dk) Op f,{a)U + J2 

i=l 

k+1 



Opf,{dk+i) + ^R'i 



^yk+l 
1=1 

This shows that 

Sn{a) = Op,,{dn) + nOnAh'-^^^) 
= Op^{dn) + On Ah'-^'^). 

In the preceding equation, the second line comes from the fact that n < log/i^^. Since 
we defined u = 1 — 26, the lemma is proved. □ 

From now on, we will always assume n < tIt for some r > fixed. We will also choose 
T small enough such that 

(3.10) r < r,,« ^ 



2r - 121oga_ 

This condition ensures in particular that Lemma 13.41 is valid, but it turns out that we 
will need a stronger condition than tF < 1/2 to complete the proof of Proposition 13. 1[ 

In order to apply the functional calculus to 5n(a), we must bound its spectrum with 
the help of Proposition 12. 1[ This is expressed in the following 
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Proposition 3.5. Define t] = — 2Tloga_ > 0. For h > small enough, we have 

Spec5„(a) C [C„,/,,2] , 

where Cn^h = e^"''°S'^- —ch}'"^^^ and c > 0. In particular, there exists a constant C > 
such that Cn,h > Ch^ , and Sn{a) has strictly positive spectrum. 

Proof. We begin by proving the result for Op;j(a„). Since we have for some ci > 

an > a^- = e2"i°S'^- > cih-'^^^"^"- , 

we will consider the symbol 6„ = a„ — e^"'°s"- > 0. Prom Lemma [331 we have 5„ € Sg+ 
and we can apply Proposition 12.11 : for any A G Spec(Op/j(a„)), there exists c > 
such that |A| > e2"'°S'^- -ch^'^^^ > cih-^^^°sa- _ ch^-25+ ^^.^g^. ^^.^^ ^ strictly 
positive spectrum, we must have — 2rloga_ < 1 — 2J^, which is satisfied if r is chosen 
according to Eq. (I3.10p . Hence, there is a constant C > such that for h small enough, 
cih~'^'^^°^°- — c/i^~^^^ > , and the lower bound is obtained. For the upper 

bound, we remark that we assumed that a+ = 1. By Proposition 12.11 it follows that any 
constant strictly bigger than 1 gives an upper bound for the spectrum. We return now 
to Sn{a)- Since we have ||5„(a) — Op;j(an)|| = 0{hy ) and v~ > rj thanks to Eq. (|3.10|) . 
we get the final result if h is small enough. □ 

Let us finish now the proof of Proposition 13. 1[ For n < n^, we begin by constructing 
a smooth function Xn,h compactly supported, equal to 1 on Spec(5ri). To do this, we 
define 

if a; < \Cn,h 



Then, the function 



Xn,h{x) = ■( 1 \ix & [C„,/i,2] 
if a; > 3 



n Xn,h{x) 
In'-X^ log X 

In 



is smooth and equal to the function 2^ log(x) on Spec5„(a), since we have shown that 
Spec(5n(a)) C [C„jj,2]. Furthermore, is compactly supported, uniformly bounded 
with respect to h. Since Cn^h ^ C/i'', the function Xn,h can easily be chosen /i^'-admissible, 
which means that In will also be /i^'-admissible. Applying the standard functional cal- 
culus, we have 

IniSnia)) = ^ logcS„(a) =^ £Sn{a) , 

where the first equality follow from the fact that £n{x) = ^logx on Spec(5„(a)). To 
complete the proof of Proposition 13.11 we compute £n(5„(a)) using both Lemma [33] and 
Corollary 12.91 For this purpose, we must check that the conditions required by this 
corollary are fulfilled. First, a„ G 5"°+, and from Eq. (IXTnll we have rj < Second, 
the remainder in Eq. (13. 7p is of order h'^ with u = 1 — 25, and we clearly have 

1 - 2b~ 
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We now apply the Corollary 12.91 to obtain 

(3.11) iSn{a) = in{Ovh{^n) + = OpJ£a„) + OuAhl , r > 0. 

Let us show that r = a~ , where a = 1 — rjTa^n- Indeed, the functional calculus states 
that 

r = min{l — 26^ — 6e, — 4e} 

where e has to be chosen in the interval min( ^~g'^^ , ^)[. Let us choose e = . Since 
1/ = 1 — 25, we have 

(3.12) r = min{l - 25+ - 6??+, - 47/+} = 1 - 25+ - 6/?+ = (1 - r/Ta^^y = a' . 
To get now 5„(a), we again apply the standard functional calculus, and we obtain 

exp {in{Sn{a))^ = exp {^^°^^r}j '= Sn{a) 

where the first equality follows again from the property of in 

on Spec(£5n). We already 

stress that this equation is essential to study the eigenvalues distribution of Mh{a,K). 
Since a„ G 5"^+, iun = 2^1oga„ € S^+. To compute 5„(a), we simply choose a smooth 
cutoff function supported in a /i— independent neighborhood of Spec(^5„(a)). The Corol- 
lary [23] is used again to get 

(3.13) Snia) = exp(£5„(a)) = OpJa„) + Ow^(/i"") . 

□ 

3.2. Radial spectral density. We begin by recalling some elementary properties of 
ergodic means with respect to the map k. Since \a\ > on T^, the function x 
log I a o k;*(x)| is continuous on for any i G N. The Birkhoff ergodic theorem then 

states that lim„^oo ^an(a^) *== (aoo{x) exists for //—almost every x. More precisely, if we 

denote Eloo(a) *== ess inf and ESoo(a) *== ess supooo for a^o = exp(^aoo), we have: 

Eloo(a) < CLoo(x) < ESoo(a) for /x — almost every x . 

In particular, if k is ergodic with respect to the Lebesgue measure /x, the Birkhoff ergodic 
theorem states that: 

(3.14) For // — a.e. X, lim loga„(a;) = / log |a| d/i = log(a) , 

and in this case, Eloo(a) = ESoo(a) = (a). 

Proof of Theorem. \1.1[ As in Proposition 12.41 we order the eigenvalues of Mh{a,K) 
and Sn{a) by decreasing moduli. Take arbitrary small e > 0, 5 > and < 7 < |5. We 

recall that 1$ = [Eloo(a) — 5, ESoo(a) + 5] and define ^n,'y- '= ^'^ \ '^n^il-y-)} so 

H{a~^{l^-)) = 1 - /i(f^„,^-) . 

Eq. (|3.14l) implies that lim„^oo Ai(^n.7-) = 0- Hence, there exists tiq € N such that 

n > no ^ /x(0 ) < 



- 6a+ -3(a) 6 - 3(a) 
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We now choose some n > uq. Applying Proposition 12.101 to the operator Sn{a), we get 
immediately: 

(3.15) E Spec(5„(a)) : s G I^} > ii{a-\l^-)) + Oh{l) ■ 

Using the Weyl inequalities, we can relate the spectrum of Sn{a) with that of Mh{a, k). 

Call dh *== #{A G Spec{Mh{a, k)) : |A| > ESoo(fl) + <5}, where eigenvalues are counted 
with their algebraic multiplicities. Hence, using Corollary 12.51 we get 

dh dh 

(3.16) dh{ESoo{a) + 6) <Y,\^k\<^Sk. 

k=l k=l 

Among the dh first (therefore, largest) eigenvalues (sj)i=i,...,d;j of 5„(a), we now distin- 
guish those which are larger than 7 + ESoo(a), and call d'^ = #{1 < i < dh, Si < 
7 + ESoo(o)} the number of remaining ones. 

Applying Proposition 12.11 to the observable a„, we are sure that for h small enough, 
Sj G Spec(5„(a)) =^ sj < 2. Hence, for h small enough, (|3.16|) induces 

(3.17) dhS < 47 + {dh - d'h){2 - ESoo(a)) ■ 
By construction, 

dh - d'h < #|s G Spec(5„(a)) : s ^ I^j , 
and using (|3.15p . we deduce 

dh-d'h<h-\f,in^^^-) + oh{i)). 

Dividing now (I3T7II by and using the preceding equation, we obtain : 

hdh < 1+ ^~^^"^ /i(l^n,7-)+Ofa(l) 

< l + l+Oh{l). 

To get the second inequality we used the condition on 7 > 0. There exists ho{e,n) such 
that for any h < ho{e,n), the last remainder is smaller than e/3. We conclude that 
Ve, 6 >0, 3/io ^ G N such that: 

(3.18) h<ho^ #{a G Spec{Mh{a, k)) : \X\ > ESoo(a) + 6^ < h'^e . 

To estimate the number of eigenvalues |A| < Eloo(a) — 6, we use the inverse propagator: 

(3.19) Mh{a,K)-' = Uh{fiy' Ophiay' = Ophia'' ° ^^)Uh{^ir' + OhmW , 

where Eq. (|2.12l) has been used for the last equality. Since Uh{K)~^ is a quantization of 
the map k~^, the right hand side has a form similar with (II. ip . with a small perturbation 
of order h. The function satisfies 

o.'^^ < < aZ^ and = (coo)^^ • 
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Hence, we also have EIoo(a^^) = (ESoo(a))~"'^, ESoo(a~^) = (EIoo(«))^"'^ • Applying the 
above results to the operator Mh{a, n)^^ , we find some /ii(e, 6) > such that 

(3.20) h < hi ^ #^^X € Spec{Mh{a, k)) : |A| < Eloo(a) - 6^ < h'^e . 
Grouping (|3.18ll and (I3.20p and taking e arbitrarily small, we obtain Eq. (II. 2|) . □ 

3.3. Angular density. From Eq. (I2.13p . we know that for h small enough, all eigenval- 
ues of Mh{a,K) satisfy |Aj| > a_/2. We can then write these eigenvalues as 

(we skip the superscript (h) for convenience). We want to show that, under the ergodicity 
assumption on the map k, the arguments 9i become homogeneously distributed over 
in the semiclassical limit. We adapt the method presented in [BDB2l[M0Kj to show that 
the same homogeneity holds for the eigenangles of the map Uh{K). The main tool is the 
study of the traces Tr(M/j(a, k)"), where n E N is taken arbitrary large but independent 
of the quantum dimension . 

Proof of Theorem \l.S\ (ii). We begin by showing a useful result concerning the trace 

Proposition 3.6. The ergodicity assumption on n implies that 

(3.21) VnGZ\{0}, lim /iTr A4(a, k)" = . 

Proof. Note here that n is fixed, independently of h. As before, we write U = Uhin) for 
convenience. Inserting products U~^U and using Egorov's property, we get: 

n-l 

(3.22) Mhia, k)"" = Ophia'jU'^ + On{h) , where a^, =^ J| a o k^^' . 

j=0 

Let eo > 0. The next step consists of exhibiting a finite open cover T'^ = Dq U Uf=i 
with the following properties : 

• Dq contains the fixed points of and has Lebesgue measure /x(-Do) < eo- Such 
a set can be constructed thanks to Proposition 12.21 

• For each i = 1, . . . , M, K^{Di)f\Di = 0. This is possible, because is continuous 
and without fixed points on \ Dq- 

Then, a partition of unity subordinated to this cover can be constructed: 

M 

i=l 

with Xi e C~(T2), Suppxi C A for < i < M. Notice that the condition on A, i > 1 
implies that XiiXi ° «;"") = 0. After quantizing this partition we write: 

M 

1=1 
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Using (I3.22P and taking the trace: 

hTT{Mhia,Kr) = /i Tr( Op Jxo) Op Ja;)C/") 

M 

4 = 1 

Since xo ^ 0, Op/j(xo)^ = ^Vh{xo) and we can perform the first trace in the right hand 
side in the basis where Op^(xo) is diagonal. This yields to 

I Tr(Op^(xo) OvhKW)\ < all Tr(Op,(xo))| + 0(1) . 
The estimates ([2JQ]l and Eq. ([22]) imply that 

/iTr(OpJxo)) = / Xodfi + O{h). 

Since n{Do) < eo, the integral on the right hand side satisfies | Jj.2 Xo dfJ-l < eo, and finally 
hTv ( OpJxo) OvMU'') = O{eo) + 0{h). 

To treat the terms i > 1, we take first e > arbitrary small and a smooth function 
Xi S such that Suppxi C Suppxi and ||xf — Xi\\c° = 0{e). Then, we use the 
symbolic calculus to write 

C/" Ovh{x^) = f/" OvniXi)^ + On Ah) + OnA^) 

= f/" Opf,{xi)U~^ Oph{x^) + On Ah) + OnJe) . 
Using the cyclicity of the trace, this gives 

/iTr(Op,(x.)Op,(a;)C/«) = /iTr([/« Op,fc)t/-"C/" Op.fc) Op,(a;)) 

+0{h) + Oie) 
= /iTr(Op,fc o Op.fc) Op, (a;)) 

+0{h) + 0{e) 
= /iTr(Op,(a; x^ ix^ o k-))[7") + + 0{s) 
= Oih) + Oie). 

In the last line we used XiiXi ° /«^") = 0. Adding up all these expressions, we finally 
obtain 

(3.23) h Tr {Mh{a, k)") = O(eo) + 0(e) + 0(/i) . 

Since this estimate holds for arbitrary small eo and e, it proves the proposition. □ 

We will now use these trace estimates, together with the information we already have on 
the radial spectral distribution (Theorem. 11.21 (i)) to prove the homogeneous distribution 
of the angles (Oi). 

Remark: This step is not obvious a priori: for a general non-normal h^^ x h^^ = 
N X N matrix M, the first few traces Tr(M") cannot, when taken alone, provide much 
information on the spectral distribution. As an example, the N x N Jordan block Jn 
of eigenvalue zero and its perturbation on the lower-left corner J7v,e = Jn + ^E^^i both 
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satisfy TV J" = for n = 1, . . . , N — 1. However, their spectra are quite different: 
Spec(J7v,£) consists in N equidistant points of modulus e^^^ . Only the trace Tr 
distinguishes these two spectra. 

By the Stone- Weierstrass theorem, any function / S C'^(S^) can be uniformly approxi- 
mated by trigonometric polynomials. Hence, for any e there exists a Fourier cutoff K £ N 
such that the truncated Fourier serie of / satisfies 

K 

(3.24) f^''He)=^he^'-'^' satisfies - < 6 . 

k=-K 

We first study the average of f^^^ over the angles {9j): 

j=l k=-K 3=1 

For each power k S [—K^K], we relate as follows the sum over the angles to the trace 

(3.25) h E e^'^^'^ = ^ Tr(M,(a, nf) + hY, ^^^^ e^'^^'^ . 

Although this relation holds for any nonsingular matrix M, it becomes useful when one 
notices that most radii Vj are close to (a) : this fact allows to show that the second term 
in the above right hand side is much smaller than the first one. 

Take 5 > arbitrary small, and denote Is *== [(a) — 5, (a) -|- 5\. From Theorem. 11.21 (i) 
we learnt that h#{rj £ Is} = 1 + Oh{l). Hence the second term in the right hand side of 
(|3.25|) can be split into: 

(3.26) h E e--^' = h E e--' + 0,(1) . 

By straightforward algebra, there exists Ck > such that 

yr€ls,yk€[-K,K], M^—I^<CkS, 

so that the left hand side in (I3.26P is bounded from above by Ck 5 + 0,(1). By summing 
over the Fourier indices k, we find 

/i-i K J. 

/i^/W(0,) = /o+ E ^/iTV(iV4(a,K)'=) + 0^j(5) + o,(l). 

1 = 1 k=-K 

Using the trace estimates of Proposition 13.61 for the traces up to |A;| = K , we thus obtain 

h E /^""H^.O = /o + Okj{5) + Oh{l) . 
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Since this is true for every S > 0, we deduce: 



We now use the estimate (|3.24l) to write: 



e being arbitrarily small, this concludes the proof. 



□ 



4. The Anosov case 



4.1. Width of the spectral distribution. If n is Anosov, one can obtain much more 
precise spectral asymptotics using dynamical information about the decay of correlations 
of classical observables under the dynamics generated by k. We will make use of prob- 
abilistic notations: a symbol a is seen as a random variable, and its value distribution 
will be denoted Pa- If we denote as before the Lebesgue measure by this distribution 
is defined for any interval / G M by : 



J I 

This is equivalent to the following property: for any continuous function / G C'(M), one 
has 





We now state a key result concerning iun when k is Anosov. 
Lemma 4.1. Set la = log(a), and 




n 



i=l 



If K is Anosov, we have 



limsupE(n3;^) < oo . 



n— >oo 




Then 
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But for K Anosov, |cjj| < e~p\^~^\ for some p > (see [Liv]). Hence 

n n 

1=1 j=i 

and the proposition follows easily. □ 

Proof of Theorem \1.3l We now have all the tools to get our main result. In this 
paragraph, we will assume that n = = £'(r log with r < Ta^n as above. It will 
be more convenient to show the following statement: for any e > and C > 0, 

(4.1) lim < i < /i^i : log|Af^| -log(a) < Cilogh-^f-^/^} = 1 . 

For h small enough, this equation is equivalent to (|1.3|) because lA^^I > a_/2. We will 
proceed exactly as in section [321 but now 6 and 7 will depend on h. 

First we define as before two positive sequences iSh)he]o,i[ and (7/i)/ig]o,i[ going to as 
/i — > 0, and such that 

and > 0. 



7/iV^log/i-i 4 

A simple choice can be made by taking 6h oc (log/i^^)^^2 and 7^ oc (log /i^^) 2^2 ^ for 
some e g]0, 

We call isi the (positive) eigenvalues of iSn{a), and define the integer dh such that 
4 = #{1 < i < h-^ : log |A,| -£a> 5h} • 
The Weyl inequalities imply 

(4.2) dh{ia + 5h)<Y,\og\\\<Y,^Si. 

i=l i=l 

Among the dh first (therefore, largest) numbers (£sj — ^a)i=i,...,d^, we now distinguish the 
d'/j first ones which are larger than 7/1, and call 

dh - d'h = #{1 < i < dh, £si -ia< 7/^} 

the number of remaining ones. Hence : 

d'h = #{!<«< h-^ : isi -£a> 7^ . 

Substracting dh^a in Eq. (|4.2p and noticing that dh — d'j^ < h^^, we get : 



d'h 



1 

(4.3) dh<^ + - yZiisi - ia) 

hOh Oh ^ 



i=l 



Recall that (.Sn — (-a = Oph{^o,n — ^a) + 0'H^{h°') for a = a > 0. From now on, a will 
denote a strictly positive constant which value may change from equation to equation. 
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Hence, the sum in the right hand side of Eq. (14. 3j) can be expressed as : 

- ia) = Trld[^,,2] (OpJ^a„ - ia) + 0{h'')) . 

The function Id[^^ 2] can easily be smoothed to give a function such that Ih{x) = 
for a; € M \ [7/j/2,3] and Ihit) = t on [7^,2]. Such a function can be clearly chosen 
u;— admissible with (log < w{h). Note also that the logarithmic decay oft;; in 
this case will always make the function Ih suitable for the functional calculus expressed 
in Proposition 12.81 and Corollary 12.91 - see the discussion at the end of ^3.11 Continuing 
from these remarks, we obtain : 

4 h' 



hdh <^ + -hTTlh{Ovh{ian-ia) + OHAh'")) 



+ ^hi:vOY>h{Th{^an - ia)) + 0{h^) 
Oh Oh 

^T + T I Mx)P.Adx) + 0{h-), 
Oh Oh Jr 

where the functional calculus with perturbations has been used. We now remark that 
Vx S SuppT/j, one has 

(4.4) Ihix)<x^V^. 

Indeed, we can clearly choose Ih such that |X^(j;)| < 1. Hence Th{x) < x, but since for 
(7 > fixed we have C(log/i~^)~2 < x for /i small enough and x G SuppX^, we get 



X < x^Vbg/Ti , 

which imply Eq. (14. 4p . Using Lemma [4m we continue from these remarks and obtain 

^/nIhix)Px„{dx) < / nx"^ P^^{dx) = E{nxl) < 00 , 

from which we conclude that 

hdh<f + 0{-=2=-) + 0{h^) ^ . 

Oh y^log h-^6h 

The theorem is completed by using the inverse map Mh{a, k)^^ , exactly as for the proof 
of Theorem [rH since || Opft(a~^ o K)C//,(a)"i - M/,(a, k)"^ || =0{h). □ 



4.2. Estimations on the number of "large" eigenvalues. In this paragraph, we 
address the question of counting the eigenvalues of Mh{a, k) outside a circle with radius 
strictly larger than (a) when k is Anosov. As we already noticed, informations on the 
eigenvalues of 5„(a) can be obtained from the function a„, via the functional calculus. 
Hence, if we want to count the eigenvalues of Mh{a,K) away from the average (a), we 
are lead to estimate the function an away from its typical value (a). More dynamically 
speaking, we are interested in large deviations results for the map k. For a G C°°(T^), 
these estimates take usually the following form |OPl IRBY] : 
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Theorem 4.2 (Large deviations). Let c > be a positive constant, and define 

ec = \og{l + c/{a)). 

If K is Anosov, there exists a function I : M.^ ]R_|_, positive, continuous and monotoni- 
cally increasing, such that 

(4.5) lim sup — log /i (x : a;„ G [£c, +oo[) < —I{ic). 

In particular, for n>l large enough, one has 

(4.6) P.„([^c,oo[) = 0(e-"^(^^)). 

We now proceed to the proof of the theorem concerning the "large" eigenvalues of 
Mh{a,t<). 

Proof of Theorem \1.4\ As before, we prove the result for iSn{a), the extension to 5„(a) 
being straightforward by taking the exponential. Define as before 

4 = #{1 <j< /i-i : log \xf^\ >£a + ic} . 

We also choose a small p > 0, and define 

d'^ = #{1 < j < h-^ : isj >£a + ec-p}. 

The Weyl inequality can be written : 

h 

dhiia + ic)<Y^ isj + {dh - d'h){ta + £c-p). 
Substracting dh{ia + do — p) on both sides, we get for p small enough 

dhP< Y^iis, - la) - d'j^iic -p)< ^{is, - la) . 
i=i i=i 

If we choose n = rir = E{Tlogh~^) and r < Ta,K as before, we can use exactly the same 
methods as above to evaluate the right hand side of the preceding equation. Let T be a 
smooth function with X = 1 on [£c — p, 2], and X = on M \ [^c — 2p, 3]. We have 

hJ2i^Sj - la) < hTvI fOp;,(^a„ - ia) + Out^ih" )) . 
i=i 

Recall that ian — la ^ Sg+ and a = l — 26 — 6r] = l — r/Ta^n- We easily check that a~ 
satisfies the condition {ii) of Corollary 12.91 since X does not depend on h. Hence, 



h p 

hy{isj-ea)< / l,,,_2p,3]{^an - ia) + 0{hn , r > 0. 

Let us show that in fact, r = a~ . By the functional calculus, we have 

r = min{l - 25+ - 6e, a~ - 4e} = min{l - 2(5+ - 6e, 1 - 2(5+ - 6r/+ - 4e} 
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where now, e > is arbitrary since X does not depend on h: this imphes immediately 
r = (j~ . Using (|4.6|) . we get 

hdh< -P.„{[ic - 2p,oo[) + 0{h''~) = 0{- e-"^(^^-2^)) + 0{h^~). 
P P 

Since p is arbitrarily small, we end up with 

(4.7) hdh = OiK"^^^^'^ +h''') = + = 0(/jmm{r/(^c-),l-r/r-,}) _ 

It is now straightforward to see that the bound is minimal if we choose 

T~ 

□ 

5. Numerical examples 

In this section, we present a numerical illustration of Theorems ll.2l and ll.3l for a simple 
example, the well known quantized cat map. 

5.1. The quantized cat map and their perturbations. We represent a point of the 
torus X G by a vector of ffi^ that we denote X = {q,p). Any matrix A G SL{2,Z) 
induces an invertible symplectic flow k : x i-^ on via a transformation of the vector 
X given by : 

X' = AX mod 1 . 

The inverse transformation is induced by A~^ , with X = A~^X' mod 1. If Tr^ > 2, 
this classical map, known as "cat map", has strong chaotic features : in particular, it has 
the Anosov property (which implies ergodicity). Any quantization Uh{A) (see [DEGj ) of 
A G SL{2,Z) with Tr^ > 2 satisfies the hypotheses required for the unitary part of the 
maps and the Egorov estimate (I2.12P turns out to be exact, i.e. it holds without 

any remainder term. 

Let us give a concrete example that will be treated numerically below. Let m G N and 
A G SL{2, Z) of the form: 

, . . _ f 2m 1 \ 

^^■^> V4m2-1 2mJ 

Then, in the position basis we have 

(5.2) Uh{A)jk = Vhexp 2mh[mk'^ - kj + mf] . 

We can also define some simple perturbations of the cat maps, by multiplying Uh{A) 
with a matrix of the form exp(^^ Op^(ff)), with H G C°°(T^) a real function, playing 
the role of a "Hamiltonian". If we denote by e^" the classical flow generated by H 

for unit time, the total classical map will be k '= e^'' oka- For reasonable choices of 
"Hamiltonians" H, k still define an Anosov maps on the torus, and the operators 

(5.3) C/^(Ai^) = e-'^OP'^(^)C/;.(^) 
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Figure 1. Spectrum of Mh{a,K) in the complex plane for = 2100. 
The dashed circles correspond to a+ and a_, while the plain circle has 
radius (a). To the left, we plot a = ai, and a = 02 to the right. 

quantize the map k. In our numerics, we have chosen m = 1, and 

OphiH) = ^sin(27rg) 

with a = 0.05. This operator is diagonal in the position representation, and for a < 0.33 
[BKj . the classical map e^'^ o is Anosov. 

Remark that because of the perturbation, Egorov property (I2.12P now holds with some 
nonzero remainder term a priori. Such perturbed cat maps do not present in general 
the numerous spectral degeneracies caracteristic for the non-perturbed cat maps |DEG| . 
hence they can be seen as a classical map with generic, strong chaotic features. 

For the damping terms, we choose two symbols of the form a{q), whose quantizations 
are then diagonal matrices with entries a{hj), j = The function ai{q) has a 

plateau ai{q) = 1 for g E [1/3,2/3], another one for q G [0,1/6] U [5/6,1] and varies 
smoothly inbetween. For the second one, we take 02(9) = 1 — ^ sin(27rq')^. Numerically, 
we have computed 

(ai)« 0.250 and (02) « 0.728. 

Fig. [U and [2] represent the spectrum of our perturbed cat map for = 2100, with 
dampings ai and 02. The spectrum stay inside an anulus delimited by Oj+ and 
stated in Eq. (I2.14p . and the clustering of the eigenvalues around {a)i is remarkable. For 
a more quantitative observation, the integrated radial and angular density of eigenvalues 
for different values of are represented in Fig. [H We check that for moduli, the curve 
jumps around (a), which denotes a maximal density around this value, and we clearly 
see the homogeneous angular repartition of the eigenvalues of Mh{a, k). 

As we can observe in Fig. [2], the width of the jumps do not depend a lot upon h, 
at least for the numerical range we have explored. This behavior could be explained by 
Theorem ll.31 which states that the speed of the clustering may be governed by log . To 
check this observation more in detail, we define the width Wh of the spectral distribution 
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Figure 2. Integrated spectral densities for the perturbed cat map. The 
radial distribution is represented to the left, the vertical bars indicate the 
value (a) for a = ai and a = a2- The angular distribution is represented 
to the right for the map Mh{ai, k). 




Figure 3. Width of the radial distribution, together with the best 2- 
parameter fits A{logh~^)~^ and the asymptotic standard errors. 



of Mh{a, k) as 



and plot Wh as a function of h- see Fig. [3l We clearly observe a decay with h~^, although 
the 2-parameter fits A(\ogh~^)~^ hints a decay slightly faster than (log Other 
numerical investigations presented in [NSj for the quantum baker map show the same 
type of decay, and a solvable quantization of the baker map allow to compute explicitely 
the width Wh, which turns to be exactly proportional to y^log/i"^. This result, together 
with the numerics presented above, seems to indicate that the bound on the decay of the 
eigenvalue distribution expressed by Theorem 11.31 may be optimal. 
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